Quench cooling under reduced gravity 



CO 

r— I 

o 



in 
<N 



o 



> 

(N 

m" 
o 

m 



D. Chatain, 1 C. Mariette, 1 V. S. Nikolayev, 1 ' 2 ^ and D. Beysens 1 ' 2 

1 Service des Basses Temperatures, UMR-E CEA / UJF-Grenoble 1, 
IN AC, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France 
2 E SEME, PMMH-ESPCI, 10, rue Vauquelin, 75231 Paris Cedex 5, France 
(Dated: March 26, 2013) 

We report quench cooling experiments performed with liquid O2 under different levels of gravity 
as provided by a magnetic gravity compensation setup. A copper disk is quenched from 270K to 
90K. It is found that when gravity is zero, the cooling time is abnormally long in comparison with 
any other gravity level. This phenomenon can be explained by the insulation effect of the gas 
surrounding the disk. The liquid subcooling is shown to drastically improve the heat exchange, thus 
reducing the cooling time (about 20 times). The effect of subcooling on the heat transfer is analyzed 
at different gravity levels. It is shown that such type of experiments cannot be used for the analysis 
of the critical heat flux (CHF) of the boiling crisis. The minimum heat flux (MHF) of boiling is 
analyzed instead. 

PACS numbers: 47.55.np, 68.03.Fg, 47.15.gm 



I. INTRODUCTION 

The rocket engines for spatial applications need to be 
able to restart under low gravity conditions to place satel- 
lites on different orbits. The solid fuel engines are cur- 
rently used for this purpose. The cryogenic fuel like hy- 
drogen and oxygen is however much more advantageous 
energetically and is currently under study for next gen- 
erations of the rocket engines. There is however an im- 
portant limitation. After some period of inactivity, the 
temperature of the fuel injectors of such engines may rise 
above the temperatures at which the fuel may be in the 
liquid state because of the solar radiation. The engine 
cannot function before the injectors are not chilled down. 
The injectors are usually cooled by putting them into 
the contact with the cryogenic fluids (i.e. with the fuel), 
i.e. by the quench cooling. One thus needs to master 
the cooling time as a function of different parameters. 
The understanding of the latter phenomenon under re- 
duced gravity conditions thus becomes important. In 
the present article we study the quenching of a hot metal 
piece in cryogenic fluids as a function of several parame- 
ters like the gravity level. 

The prediction of cooling time necessitates the knowl- 
edge of the heat transfer coefficient which depends on 
temperature. To measure the influence of the grav- 
ity level on the value of the heat transfer coefficient, 
we have performed experiments in a liquid oxygen-filled 
cell. A copper disc was quickly immersed into this 
bath and cooled down (quenching method). The Earth 
gravity was compensated by the magnetic force in the 
OLGA (Oxygen Low Gravity Apparatus) facility located 
at CEA/Grenoble in France. 

The objective of this article is three-fold. We aim to 
show that (i) one cannot measure the critical heat flux 



of the boiling crisis in a quenching experiment and (ii) 
the subcooling improves greatly the cooling efficiency in 
microgravity that is otherwise very poor. Finally, we 
measure the dependence of the minimum flux of boiling 
(i. e. the threshold of the departure from film to nucleate 
boiling) as a function of gravity. 



A. The quench cooling 
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FIG. 1. Temperature evolution of a copper disk initially at 
about 270K dipped into O2 bath at 90K under Earth gravity. 

Consider what happens when a hot sample is dipped 
into a liquid. Fig. [T] presents a typical temperature evolu- 
tion of a copper disk initially at about 270K dipped into 
a liquid O2 bath at 90K under Earth gravity. 

At a first stage (a-b zone), the temperature decreases 
slowly because the solid surface is completely covered by 
a layer of vapor which insulates the solid from the liq- 
uid. This is the film boiling regime, during which the 
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vapor bubbles are generated at the vapor-liquid bound- 
ary. When the surface temperature becomes low enough, 
the liquid may come into direct contact with the sample 
and the bubbles start nucleating directly at the solid sur- 
face. This signifies the beginning of the nucleate boiling 
regime which provides much more efficient heat transfer 
and the temperature decreases very quickly (b-c zone). 
At the final stage, the solid temperature is not any more 
sufficiently high to induce boiling; the solid temperature 
approaches that of the liquid O2 bath due to convection. 
The evolution of heat transfer regimes can be traced on 
the Nukiyama diagram (Fig. [2]) presenting the heat flux 
q from the heater as a function of the sample tempera- 
ture T. The system evolves from the right (film boiling 
regime) to the left in Fig. O The moment where the di- 
rect solid-liquid contact begins to occur corresponds to 
the minimum heat flux (MHF). From this moment, the 
transition from film to nucleate boiling starts. The line 
along which the system evolves depends on the thermal 
inertia of solid. If it was vanishingly small, a transition to 
the nucleate boiling regime would occur instantaneously 
(leftwise directed arrow in Fig. [2]). In reality, a finite time 
is necessary to cool down the solid and the system evolves 
along the dotted curve. The heat flux attains a maximum 
value at a point where no part of the solid is covered by 
the persisting vapor film. The nucleate boiling regime 
occurs from that moment on. It is very important to 
distinguish the quenching curve from the boiling (solid) 
curve. 
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FIG. 2. (Color online) Schematic heat transfer curves for 
stationary boiling experiments and transient quenching. 



ies for the metallurgy (steel quenching), it became clear 
that only MHF may be obtained by this method [J], the 
maximum flux value depending strongly on the sample 
thermal inertia. The actual CHF may only be obtained 
for the samples of the large enough thermal inertia. How- 
ever the CHF may hardly be studied in that case, which 
is evident from the following consideration. 

During the cooling, the heat Q transferred from the 
sample to the liquid bath during the time At may be 
expressed as 

Q = mC P (T)(T(t)-T(t + At)), (1) 

where m is the mass of the disk, T its temperature, and 
t is the time. The specific heat of copper Cp depends on 

C P (T) = -1.355 x 10" 7 T 4 + 1.303 x 10" 4 T 3 
-4.798 x 10" 2 T 2 + 8.331T - 217.964 

The transferred heat may also be expressed with the heat 
transfer coefficient ft, 

Q = hS (T — T L ) At, (2) 

where S is the total surface of the sample and Tl is the 
temperature of the liquid bath. By equating these ex- 
pressions, one finds 

= mC P [T(t)][T(t)-T(t + At)} 
SAt[T(t) - T L ] 

and the heat flux q = Q/ (SAt). 

This method of measurements is valid only if the tem- 
perature of the sample is uniform during the cooling. For 
this, the Biot number 



must be lower than 0.1 [5]. Here L c is the characteristic 
length (a half thickness of the disk); A is the copper ther- 
mal conductivity. It has been checked that the criterion 
Bi < 0.1 holds for all described experiments. 

One can see now that the requirements of small Bi < 
0.1 and large thermal inertia of the sample are contra- 
dictory. This explains why the quench method is hardly 
suitable for the CHF study. The objective of this article 
is the analysis of the MHF rather than of CHF. Besides 
we study the cooling dynamics and the impact of sub- 
cooling on it at various gravity levels. 



The maximum heat flux (called critical heat flux, CHF) 
can only be recovered in the stationary boiling experi- 
ments (see, e.g. [l|) performed at fixed heat flux during 
each experiment. Till 1960s, the quenching experiments 
were considered to be suitable for the CHF determina- 
tion [2, 3]. With the development of the quenching stud- 



II. EXPERIMENTAL APPARATUS 

Experiments were performed in the OLGA facility 
(Figl3]). It is equipped with a superconductive coil and 
creates a magnetic force strong enough to compensate 
the buoyancy force in O2. 
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FIG. 3. (Color online) A sketch and a general view of the 
OLGA facility. 




FIG. 4. (Color online) A sketch of the cell and a photo of its 
lower part (the upper part of the cell is removed). The disk 
is shown in the upper position on the sketch and in the lower 
on the photo. 



Oxygen is contained in the experimental cell (FigHJ) 
placed in an anti-cryostat inserted into the 300 mm bore 
of the superconductive coil. Two endoscopes are used for 
observation: one for the light and another for the camera. 

The lower part of the experimental cell is made of a 
sapphire tube of 30 mm diameter and 100 mm length 
closed by two copper flanges equipped with heat exchang- 
ers, inside which the gaseous cold helium may flow. The 
temperature of the two flanges is controlled within 0.0 IK 
with heaters. The upper part of the cell is the stainless 
steel tube. The pressure in the cell is measured with 
an accuracy of 10 mb. Liquid oxygen is obtained by 
condensing pure O2 gas (99.995%) in the cell. The cell 
is brought to the desired temperature by regulating the 
power of heaters and helium flow rate in the heat ex- 



changers that are situated at the top and at the bottom 
of the transparent part of the cell. 

The disk is made of pure copper. Its diameter is 
20.02 mm and its thickness is 2.98 mm. The weight 
of the disk is 7.68 g, its surface in contact with liquid 
is S = 8.07 cm 2 and its RMS roughness is 0.8 jam. A 
hole was drilled to introduce a platinum thermometer 
(supplier: Heinz Messtechnic, model W60-50, class A; 
diameter 1.8 mm, length 11 mm) in its heart. After each 
quench experiment, the disk may be risen into a hot zone 
that situates 0.6 m above the transparent part of the cell 
with a screw/nut system driven by a stepper motor. The 
hot zone of the cell is regulated at 270 K. Once the disk 
attains this temperature, the grip is loosened again and 
the disk drops into the liquid. Meanwhile, the disk tem- 
perature is acquired at a frequency of 100 Hz. 

III. GRAVITY COMPENSATION 

The principle of magnetic compensation of gravity in 
a pure non-magnetic substance is an application of the 
magnetic volume force opposite to the Earth gravity 

This force per unit volume is given by the expression 

Fra = |^V(B 2 ), (3) 

where B is the magnetic field, fio is the magnetic perme- 
ability of vacuum and Xm the magnetic susceptibility of 
the material. 

Since Xm is proportional to the density p of the non- 
magnetic substance, one may introduce a quantity a = 
Xm/ P independent of the density. It means that the buoy- 
ancy force is compensated if one applies the field with 
dB 2 jdz = 2jiQg/a 1 where g is the Earth gravity acceler- 
ation and z is the vertical coordinate. For O2 at 90K, 
2fi g/a ~ 8 T 2 /m. These conditions are achieved in 
OLGA facility near the bottom of the superconductive 
coil. 

It was demonstrated [9[ that the magnetic force ((3]) 
cannot be made constant in a volume. In our experimen- 
tal cell, Earth gravity may be exactly compensated at a 
single point called complete compensation point. Such 
conditions are called hereafter Og and the corresponding 
coil current is denoted I$ g . In the vicinity of this point, 
a residual acceleration acts on O2 [10]. As the magnetic 
force is a way too small to compensate the gravity for cop- 
per, the disk falls with the Earth gravity acceleration. Its 
stop point situates at z = —242 mm with respect to the 
center of the coil. At this location, Fig. Ogives the calcu- 
lated iso- values of residual gravity acceleration modulus 
at Og conditions. One can see that at the edge of the 
disk, the residual acceleration is smaller than 2% of g. 

To achieve Og at the required place inside the cell, a 
single bubble is generated by the boiling process with 
a heater situated at the bottom of the cell. Then the 
current in the coil is adjusted until the bubble center 
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FIG. 5. (Color online) Iso lines of the residual acceleration 
modulus in % of g inside the cell of OLGA at the ending point 
of the disk fall at Og. 



coincides with the optimal previously calculated point 
(Fig. ED. 

The effective gravity level |8[ may be varied by chang- 
ing the coil current /. It is known [11| that the effective 
gravity level g* at the geometrical point of complete com- 
pensation is given by the equation 




FIG. 6. Image of O2 bubble at Og inside the cell at equilibrium 
as seen in the periscope mirror used to adjust the position of 
the exact compensation point with respect to the cell. For 
this case, the temperature is 92.24K and the current in the 
coil is Io g = 241 A. The elongated shape of the bubble is due 
to the weak coupling of the magnetic field with the gas-liquid 
interface shape [8|- 

The magnetic susceptibility of oxygen varies with tem- 
perature, and Io g depends weakly on the liquid temper- 
ature. For the isothermal conditions, the currents for 
the complete gravity compensation are Iq 9 = 239.7 A at 
90.07K (pressure of 1 bar, at saturation) and Io g = 241 A 
at 92.24K (subcooling 5K for the pressure of 2 bars). 
When the disk is dropped into liquid O2, the liquid is 
slightly heated up, resulting in a small decrease of the 
magnetic susceptibility. By observing the movement of 
the gas generated by boiling around the disk (Fig. [TJ), we 
could see that the current had to be slightly increased to 
obtain Og in the cell during the disk cooling. In this case, 
log = 242.1 A at 90.07K, and I 0g = 245.18 A at 92.24K. 
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FIG. 7. (Color online) Film boiling at different gravity (its 
direction is shown with vertical arrows) at initial liquid O2 
temperature (before disk drop) 92.24K. a) I = 241 A: the 
residual gravity is directed downwards and the gas surround- 
ing the disk slowly rises, b) / = 245.18 A: the gas recondenses 
without rising and surrounds the disk during all the cooling 
time (zero gravity), c) / = 252.7 A: the residual gravity is 
directed upwards and the gas goes downwards. 



IV. EXPERIMENTAL RESULTS 
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FIG. 8. (Color online) Evolution of the copper disk temper- 
ature dipped into O2 bath at different gravity levels for the 
saturated case. 



Quenching experiments were performed for two O2 
thermodynamic conditions: saturated O2 and subcooled 
O2 with respect to its saturation temperature for given 
pressure that was controlled independently since the cell 
was connected to a large oxygen-filled vessel. 



A. Case of saturated liquid 

For saturated O2 experiments, the temperature of the 
cell was regulated at 90.07K and the vessel pressure was 
maintained at 1 bar (which is the saturation pressure 
corresponding to this temperature). 15 quenches were 
performed. 
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1. Cooling curves of the disk 

Fig. [8] shows the cooling curves of the disk for different 
levels of gravity. We can see that the lower the gravity, 
the larger is the cooling time. The Og cooling time is very 
large with respect to other cases and the corresponding 
curve is truncated: it took more than 10 minutes to cool 
down the disk. The complete curve is shown in the inset 
to Fig. [T3l below. 



2. Heat transfer curves 



Using the above described method, we have calcu- 
lated the heat transfer curves for different levels of grav- 
ity (Fig. [9]). The film boiling points were calculated 
by making a centered moving average on 50 points for 
the film boiling and from 10 to 20 points for the transi- 
tional/nucleate boiling. One can clearly see the influence 
of the gravity level on the transferred heat flux in film 
boiling regime and on the maximum and minimum heat 
flux values. 
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FIG. 9. (Color online) Heat transfer curves for saturated O2 
quench cooling at 1 bar and different gravity levels. 



Fig. [10] is a comparison of Kutateladze and Breen & 
Westwater correlations [12| with our experiment under 
lg. One can see that the quench method gives the max- 
imal heat flux ~ 12 W/cm 2 ), which is two times lower 
than the CHF value 24.9 W/cm 2 obtained with the Ku- 
tateladze correlation. The latter value reasonably agrees 
with the experimental value (22-23 W/cm 2 ) obtained in 
pool boiling experiments [l3| . The low value of the maxi- 
mum heat flux reflects the transient nature of the quench 
cooling as discussed in Fig. [2f sec. II A|) . 
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FIG. 10. (Color online) Comparison of the Kutateladze (bro- 
ken line) and Breen& Westwater (circles) correlations with our 
data from Fig. [9] corresponding to lg (solid line). 
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FIG. 11. (Color online) Sketch showing the principle used for 
subcooling experiments. 



B. Case of subcooled liquid 

To measure the subcooling influence, experiments 
where performed for several subcooling values. The 
subcooling was obtained with the following method (cf. 
Fig. E]) . The temperatures of two lower heat exchangers 
of the cell (cf. Fig. H]) are regulated at the same value 
T\. The temperature T2 of the upper heat exchanger is 
regulated at a higher value equal to the saturation tem- 
perature corresponding to the gas pressure. Under these 
conditions, the equilibrium level of liquid forms in the 
region of cell maintained at the temperature T2. A ther- 
mal gradient exists in the liquid between two upper heat 
exchangers. The liquid is isothermal at T\ in between 
two lower heat exchanges. The final drop point of disk is 
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situated in this initially isothermal region. 

To achieve the 5K subcooling, T x = 92.24K and 
T2 = 97.24K. The latter temperature corresponds to the 
pressure of 2 bars. For the 7K subcooling (see below), 
Ti = 90.24K. 



1. Cooling curves 

Fig. [12] shows the raw cooling curves of the disk for 6 
levels of gravity. The temperature of cell was regulated 
at 92.24K and the pressure was adjusted at 2 bars which 
results in a subcooling of 5K. The temperature acquisi- 
tion frequency was 100 Hz. Similarly to the saturated 
case, the lower the gravity level, the longer is the cooling 
time. Unlike the saturated case, the Og cooling time is 
only twice larger than that of the lg case. 
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FIG. 12. (Color online) Evolution of the disk temperature at 
different gravity levels for the 5K sub cooled case. 



2. Impact of subcooling 

Fig. [13] shows that the subcooling reduces drastically 
(20 times) the general cooling time at Og while the ef- 
fect at any nonzero gravity is much weaker (cf. Figs. [8] 
and [T2]). Such a reduction may be explained as follows. 
Fig. 03 shows the formation at Og of a vapor bubble that 
surrounds the solid and completely insulates it from the 
liquid. The bubble can recondense if the liquid is sub- 
cooled and cannot recondense if the liquid is saturated. 
For the saturated case, the heat exchange is thus pro- 
duced via the weak heat conduction through the vapor 
and through the almost non-conductive at cryogenic tem- 
peratures stainless steel rod that supports the disk (cf. 
Fig. 3]). When the liquid is subcooled, the recondensa- 
tion of the vapor at the external bubble interface is very 
efficient and the gas layer surrounding the disk is much 



thinner. This effect improves the heat exchange thus fa- 
voring cooling. 
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FIG. 13. (Color online) Evolution of the disk temperature at 
different subcooling levels at Og. The same curves are shown 
in the inset in the semi logarithmic scale. 



3. Heat transfer curves 

We have plotted the heat transfer curves for 6 levels 
of gravity (Fig. fl4|) using the same method as the one 
used for the saturated case. One can see that unlike the 
saturated quench, the maximum and the minimum heat 
flux values are both nonzero for the Og case. 
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FIG. 14. (Color online) Heat transfer curves for the 5K sub- 
cooling case at different gravity levels. 
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FIG. 15. (Color online) Minimum heat flux versus gravity 
level for saturated and subcooled 02- Inset: same in log- log 
scale; the 1/4 slope is shown for comparison. 



C. Minimum heat flux 



Fig. [15] shows the minimum heat flux as a function 
of the gravity level. The minimum value was obtained 
by fitting the raw data with a second degree polyno- 
mial. The MHF occurrence is usually associated with 
the Taylor instability of the vapor film [14], ll5| so that 
Qmhf ~ g 1 ^ should be satisfied. Fig. [T5l shows an expo- 
nent close to this value. 



V. CONCLUSION 

Quench cooling has been analyzed under different grav- 
ity levels. Earth gravity was compensated with magnetic 
forces. Measurements were performed by quenching a 
copper disk of 20 mm diameter and 3 mm thickness from 
270K to 90K. The analysis of different boiling regimes 
has shown that the study of the critical heat flux (CHF) 
via quenching experiments is hardly possible because of 
the transient nature of the quench cooling: the maximum 
heat flux value is twice smaller than the CHF measured 
by other authors in the stationary boiling experiments. 

It has been shown that cooling via quenching in micro- 
gravity takes usually exceedingly long time because the 
vapor generated during the boiling envelopes completely 
the cooled body (film boiling regime), thus thermally in- 
sulating it from the liquid. It is shown that the artificial 
subcooling (for instance, by transient pressurization) can 
drastically accelerate the cooling. This speeding up is ex- 
plained by the recondensation of the vapor envelope at 
the vapor-liquid interface in case of subcooling. 

The dependence of the heat transfer and, in particular, 
of the minimum heat flux of boiling (MHF) on gravity is 
analyzed. As expected, the heat transfer improves with 
the gravity level. The MHF gravity dependence shows 
an exponent close to that of the Zuber theory (exponent 
1/4). 
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